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INTRODUCTION 


In  an  earlier  paper,  Klett1  presented  a  stable  solution  to  the  single¬ 
scattering  lidar  equation 


S(r)  =  ln[p(r)r2]  =  ln[c  C  ]  +  k  lnc(r)  -  2  /  o(r)dr. 


,  CTAr 

In  this  equation  P(r)  is  the  power  received  at  a  range  r  and  Cj  is  p^ — — 

where  PQ  is  the  transmitted  power,  c  is  the  speed  of  light,  x  is  the  pulse 

width,  A  is  the  receiver  area  and  F  is  the  overall  system  response.  The 

constants  C2  and  k  are  atmospheric  dependent  parameters  which  relate  the 

coefficients  of  backscatter,  8  (r),  and  extinction,  a(r) ,  according  to  the 

it 

power  law  relationship1 


P*(r)  “  c2°(r) 


From  Equation  (1)  a  differential  equation  (independent  of  C1  and  C2 ) 


<*SU)  «  *  dCT(r)  _  2o(r) 

dr  d(r)  dr 


can  be  obtained  which  has  the  well  known  solution 


o(r)  - 


exp{[s(r)  -  S(r  )]/k} 


d(rQ)  1  /  exp{[s<r)  -  S(rQ)]/k}dr 

ro 


where  a(r  )  is  the  extinction  coefficient  at  a  range  rQ.  This  range  is 
normally  chosen  where  the  transmitted  beam  and  receiver  field-of-view 
overlap.  In  this  equation  the  extinction  is  detemined  by  the  ratio  of  two 


i 


numbers  which  each  become  smaller  as  r  increases.  In  particular,  for  good 

visibilities  the  denominator  becomes  the  difference  between  large  numbers  and 

small  flucuations  in  S(r)  can  cause  large  oscillations  in  o( r) .  Klett1  also 

discusses  the  instabilities  in  <J( r)  caused  by  uncertainties  in  estimating 

o(rQ) .  He  proposed  a  more  stable  solution  to  Equation  (3)  for  r<r^  rather 

than  r >r  .  This  new  solution  is 
o 

exp{[s<r)  -  S(r  ,) ]/k} 

c(r) - - -  (5) 

r 

o(rf)  1  +  /  fexp{[s(r)  -  S(rf )  ]/k}dr 

r 

where  cr(r^)  is  the  extinction  coefficient  at  range  r^.  In  this  case  it  is  not 
possible  for  the  denominator  to  become  zero  or  negative  and  0(r)  is 
determined  by  the  ratio  of  two  numbers  which  become  larger  as  r  decreases  from 

rf 

Recently,  Ferguson  and  Stephens2  described  a  novel  approach  in  deter¬ 
mining  o(r)  using  both  Equations  (4)  and  (5).  In  their  algorithm  an  iterative 
scheme  is  used  to  select  the  appropriate  value  of  0(r)  using  Equation  (5) .  A 
value  of  o(r)  is  determined  such  that  the  value  of  Sc(r)  calculated  from 
Equation  (1)  matches  the  measured  value  Sffl(r).  The  iteration  procedure  is 
initiated  at  a  close-in  range  where  the  returned  signal  is  well  above  the 
system  noise  and  requires  a  minimum  amount  of  computational  time.  To 
demonstrate  the  utility  of  their  algorithm,  Ferguson  and  Stephens2  presented 
examples  of  calculated  and  measured  values  of  S(r)  normalized  to  the  product 
0^2*  RMS  differences  of  less  than  10”4  over  the  range  rQ  to  rf  indicate  the 

value  of  a(r  )  to  be  very  accurate.  This  value  was  then  used  as  <j(r  )  in 
f  o 

Equation  (4)  to  determine  o(r)  for  r>r^.  In  this  case  it  is  possible  to 
extend  c(r)  out  to  ranges  where  S,^*)  disappears  into  system  noise.  However, 


it  Bast  be  realised  that  the  solution  for  o(r)  is  not  unique  unless  the  values 
of  C2  and  k  along  the  propagation  path  are  known.  The  proper  choice  of  these 
parameters  is  a  critical  problem  in  interpretation  of  lidar  returns.  In  his 
original  paper,  Klett1  showed  Equation  (5)  to  be  rather  insensitive  to  changes 
in  k  for  highly  turbid  atmospheres.  Little  attention  has  been  given  to  the 
sensitivity  of  the  inversion  algorithms  to  changes  in  C2  which  is  usually 
assumed  to  be  invariant.  However,  from  the  work  of  Barteneva3,  a  change 
greater  than  an  order  of  magnitude  can  be  inferred  in  the  value  of  C2  between 
clear  air  and  fog  conditions.  Fitzgerald^  has  also  pointed  out  that  a  power- 
law  relationship  between  backscatter  and  extinction  is  only  valid  for  relative 
humidities  greater  than  about  80  percent  and,  even  then,  it  is  dependent  upon 
the  air  mass  characteristics. 

In  this  paper,  we  investigate  the  sensitivity  of  the  inversion  algorithm 
to  uncertainties  in  C2  and  k.  For  this  study  we  utilized  data  acquired  in 
September  1982  at  a  research  site  located  near  the  Baltic  Sea.  The  data  were 
obtained  using  a  handheld  battery  operated  lidar,  termed  a  "visioceilometer ," 
developed  by  the  U.S.  Army  for  measuring  atmospheric  visibility  and  cloud 
heights.5  The  optical  unit  is  a  modified  AN/GVS-5  laser  rangefinder  which 
emits  a  10-m  J,  6-nsec  pulse  at  1.06^m.  The  signal  processing  unit  clocks  the 
return  signal  through  a  transient  recorder  at  a  20-MHz  rate  giving  a  7.5  m 
sampling  interval.  The  digitized  results  are  transferred  to  a  microprocessor 
and  then  to  a  tape  recorder  for  off-line  processing  by  a  Hewlett-Packard  9845B 
computer. 


INVERSION  EXAMPLES  AND  ANALYSIS 

For  this  study  we  have  chosen  a  lidar  return  obtained  during  reduced 
visibility  conditions  for  which  the  Klett  inversion  (Equation  5)  rapidly 
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converges.  Figure  1  is  a  plot  of  S(r)  versus  range  determined  from  a  lidar 
return.  The  initial  slope  of  the  curve  shows  a  nearly  monotonic  decrease  with 
range  and  indicates  a  homogeneous  atmosphere  out  to  nearly  500  m.  In  this 
case,  (from  Equation  1)  the  extinction  coefficient  is  given  by  -  • 

Beyond  500  m  the  atmosphere  appears  less  homogeneous  and  may  represent  a 
variation  with  altitude  as  the  lidar  was  elevated  3° .  Here  we  will  limit  the 
analysis  to  the  homogeneous  portion  of  the  return  between  ranges  rQ  “  112.5  m 
and  »  412.5  m.  For  this  portion  of  the  curve,  a  least  squares  straight 
line  fit  (with  a  correlation  coefficient  of  0.956)  gives 

S(r)  -  -1.68r  -  4.98  (6) 

from  which  the  homogeneous  extinction  coefficient  Oq  is  determined  to  be  0.84 
km-1.  By  equating  Equation  (6)  to  Equation  (1)  for  the  homogeneous  case  we 
find  (for  k  «  1) 

ln^C  ]  -  -4.8.  (7) 

If  the  system  constant  C,  is  known,  then  C2  is  determined  for  this  special 
case.  However,  since  the  overall  system  response  of  the  visioceilometer  is 
not  known,  we  will  in  fact  determine  the  sensitivity  of  the  algorithm  to 

uncertainties  in  the  parameter  In [c^C  ]  which  are  directly  related  to 

variations  in  C2  (for  fixed  values  of  k)  by  the  expression 

lnf^c’]  -  ln[c’/C2]  -  4.8  (8) 
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S(r)  =  -1 .68r  -  4.98 
(p  =  .956) 


RANGE  (km) 

Figure  1.  S(r)  versus  range  determined  from 

a  lidar  return. 
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I  • 

where  differs  from  the  correct  value*  C2#  In  Figure  2  the  value  of  C 2  is 
allowed  to  change  by  ±40%  from  C2  with  resulting  changes  in  lnfc^C^]  from  +7% 
to  -10%. 

Figures  3A  to  I  are  plots  of  the  extinction  coefficients  calculated  using 
the  values  of  ln[ciC2]  shown  in  each  figure  and  the  S(r)  values  from  Figure  1 
with  k  -  1.  For  these  calculations  the  RMS  differences  between  the  measured 
and  calculated  values  of  S(r)  ranged  between  10-2  and  10“4.  The  extinction 
coefficient  in  Figure  3E  (ln[c,jC2]  =  -4.79)  is  nearly  constant  out  to  500  m 
and  is  in  close  agreement  with  that  determined  from  the  slope  of  S(r)  in 

I 

Figure  1.  However,  for  the  cases  where  c2/c2  <  1  (Figures  3A  to  D) ,  the 
calculated  extinctions  increase  with  range  to  a  value  near  4  km-1  in  Figure 

1 

3A.  When  c  /C  >1  (Figures  3F  to  I)  the  extinction  coefficients  decrease 

2  2 

with  range  and  tend  towards  zero  in  Figure  31.  These  tendencies  are  similiar 
to  the  singularities  and  zeros  discussed  by  Klett1  when  the  inversion  is 
started  at  the  beginning  of  the  return  signal. 

Also  listed  in  the  figures  are  the  visibilities  determined  from  the 
relation 


VIS 


3.912  (r  -  r  ) 
f  o 


/  fc(r)dr 
r 

o 


(9) 


In  Figure  4  the  calculated  visibilities  are  plotted  for  the  relative  changes 
in  c2  and  differing  values  of  k.  In  general,  these  calculations  show  the 
visibility  to  be  rather  insensitive  to  k.  However,  we  have  found  the 


extinction  coefficients  to  be  sensitive  to  small  changes  in  k  when  c(r)  is  at 
or  approaching  a  singularity  as  in  Figures  3A  and  3B.  The  sensitivity  of 


rigor*  2.  Changes  in  the  paraneter  lnfc^C^]  with  changes 
in  the  value  of  Cj  appropriate  for  Equation  7. 


7 


LN(C,Cj)  =  -5 
VIS  =  2.28  KM 


VISIBILITY,  km 


visibility  to  relative  uncertainties  in  C2  is  seen  in  Figure  4  where  a  20 
percent  uncertainty  can  cause  nearly  a  30  percent  change  in  visibility. 


The  variations  of  extinction  coefficie;  '•s  with  range  can  be  explained  for 
a  uniform  atmosphere  [a(r)  “  Oq]  if  we  use  values  of  and  o' (r)  to  evaluate 
the  closeness  of  the  measured  and  calculated  values  of  S(r),  i.e.. 


Sm<r> 


SC(D 


(10) 


or 


ln[ciC2]  +  k  ln(Oo)  -  2  CQr  =  lnfc,^]  +  k  lno*(r) 


-  2  a’  (r  )  r  -2  /  a'  (r)dr. 
o  o 

r 

o 


Rearranging  gives 


(11) 


-1  n[c2/C2]  +  k  ln[an/<j'(r)]  -  2 o^r  +  2  o'  (rj  r 


o  o 


+2  /  o' (r)dr  =  0 
r 


For  a  uniform  atmosphere, 


(12) 


S  (r)  -  S  (r  )  -  -2a  (r  -  r  ) 
at  mo  o  o 


(13) 


and  o' (r)  are  then  determined  from  Equation  (4)  to  be  (for  k  »  1) 


a'  (r) 


o'  (r  ) 
o 


exp[-2a0(r  -  rQ) ] 


r 

-  2 J  exp[-2oQ(r 
r 

o 


(14) 
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o 


2 


For  X  >  1  we  find  that  <?'  (r)+  0  as  r+  “.  Denoting  the  left  and  right  sides  of 

Equation  (17)  as  f 1 ( X)  and  f^ ( X) »  respectively,  a  graphical  solution  to  the 

equation  is  presented  in  Figure  5  for  o  =  0.85  km  From  the  curves  it  is 

o 

•  i 

seen  that  x  >  1  when  C2/C2  >  1  and  X  <  1  for  C2/C2  <1.  An  example  of  a 
numerical  solution  of  Equation  (17)  is  shown  in  Figure  6  which  illustrates  the 
behavior  of  o' (r)  with  range  for  uncertainties  in  C2>  if  c2  is  under¬ 
estimated,  o’  ( r)  tends  to  increase  without  bound.  If  c2  is  underestimated  by 
20  percent,  o'  (r)  differs  from  the  known  oq  by  a  factor  of  4  at  a  range  less 
than  1  km.  If  C2  is  over  estimated,  fl'(r)  tends  to  zero.  These  trends  in 
o’  (r)  are  identical  to  those  illustrated  from  the  measured  data  in  Figure  3. 


(X)* 


•<J>  o 


o' - 1 - 1 - 1 - JL 

0  1  2 


RANGE,  km 

Figure  6.  Numerical  solution  of  Equation  (17)  for  a’  (r)  versus 
range  for  relative  changes  in  C2  • 


14 


CONCLUSIONS 


This  analysis  has  determined  that  range  dependent  extinction  coefficients 
inferred  from  lidar  returns  are  extremely  sensitive  to  the  value  of  the  con¬ 
stant  relating  backscatter  and  extinction  for  a  reduced  visibility  condition* 
The  instabilities  will  also  occur  for  better  visibilities  but  at  longer 
ranges*  A  correct  interpretation  of  atmospheric  structure  deduced  from  lidar 
data  requires  precise  knowledge  of  C2  along  the  path.  It  is  also  important  to 
realize  that  changes  in  extinction  coefficients  occur  with  changes  in  the 
product  C^c2*  Uncertainties  or  instabilities  in  the  system  constant  C^  can 
result  in  similar  instabilities  in  the  extinction  coefficient. 
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